Robust-to-loss entanglement generation using a quantum 
plasmonic nanoparticle array 



Changhyoup Lee, 1 2 Mark Tame, 3 Changsuk Noh, 1 James Lim, 2 4 5 
Stefan A. Maier, 3 Jinhyoung Lee, 2 4 and Dimitris G. Angelakis 1 6 

Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, 
Singapore 117543 

2 Department of Physics, Hanyang University, Seoul 133-791, Korea 

3 EXSS, The Blackett Laboratory, Imperial College London, Prince Consort Road, SW7 

2BW, United Kingdom 

4 Center for Macroscopic Quantum Control, Seoul National University, Seoul 151-742, Korea 
5 Research Institute for Natural Sciences, Hanyang University, Seoul 133-791, Korea 
6 Science Department, Technical University of Crete, Chania, Crete, Greece, 73100 

E-mail: chang doll i@gmail . com and dimitris.angelakis@gmail.com 

Abstract. We introduce a scheme for generating entanglement between two quantum dots 
using a plasmonic waveguide made from an array of metal nanoparticles. We show that the 
scheme is robust to loss, enabling it to work over long distance plasmonic nanoparticle arrays, 
as well as in the presence of other imperfections such as the detuning of the energy levels of 
the quantum dots. The scheme represents an alternative strategy to the previously introduced 
dissipative driven schemes for generating entanglement in plasmonic systems. Here, the 
entanglement is generated by using dipole-induced interference effects and detection-based 
postselection. Thus, contrary to the widely held view that loss is major problem for quantum 
plasmonic systems, we provide a robust-to-loss entanglement generation scheme that could be 
used as a versatile building block for quantum state engineering and control at the nanoscale. 
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1. Introduction 

Quantum plasmonics is a rapidly emerging field that offers new opportunities for investigating 
quantum optics at the nanoscale JD El 0. Here, novel capabilities in the way the 
electromagnetic field can be localised 0) and manipulated [5] open up the prospect of 
miniaturisation, scalability and strong coherent coupling with single-emitter systems J61, 
beyond the limits of conventional photonic systems 0. In particular, with the advancement 
of nanofabrication and characterisation technologies, metal nanoparticles have been attracting 
considerable attention as they allow a flexible approach to reaching a high confinement of 
optical fields [[8l |9l [10l fTTTl . and it has recently been suggested to use them for building 
compact on-chip quantum plasmonic networks operating at the nanoscale [Q21 \T3\ fT4ll . 
However, Ohmic loss in the metals that support plasmonic excitations is a major obstacle for 
realising plasmonic quantum information processing and quantum control lfT5ll . This energy 
dissipation process induces decoherence in the system and limits the performance of a given 
task. For example, it has been shown that the optimal distance for quantum state transfer 
in an array of nanoparticles is severely limited. lfT2ll . In order to overcome the problem 
of loss in plasmonic systems researchers have begun to consider employing various types 
of gain media lfT6l or metamaterials [17], while trying to keep the high field confinement 
characteristic offered by plasmonic systems. A different approach has been to tailor the system 
dynamics such that the steady state of the system yields useful states, the so-called driven- 
dissipative approach H8l[T9]]. Despite this important progress in techniques to overcome loss 
in plasmonic systems, implementing quantum information processing based on an array of 
metal nanoparticles remains elusive. 

In this paper, as a first step towards building quantum plasmonic networks based on metal 
nanoparticles, we propose an alternative and practical scheme for generating entanglement, 
between two distant quantum dots (QDs), using the plasmonic modes of a metal nanoparticle 
array. We show that one can achieve entangled QD states with high fidelities even in the 
presence of large losses from the metal nanoparticles and inhomogeneous broadenings of the 
QDs. Robustness against metal loss of the entanglement (or equivalently against an increase in 
the length of the nanoparticle array) alleviates the length limit of the array, opening up its use 
for compact nanoscale quantum networks, and for providing opportunities in the construction 
of quantum plasmonic devices on scales far below the diffraction limit. The robustness of the 
entanglement against detunings of the QDs allows the scheme to work even when the QDs 
have different resonant frequencies, where the operating frequency to be detected is far off- 
resonant with the QDs transitions. This is an important practical consideration. It ultimately 
stems from the fact that only classical lights are required to generate the entanglement between 
two distant QDs. The scheme we introduce is based on a quantum interference phenomena 
introduced by Waks and Vuckovic [20 j called dipole-induced transparency (DIT). We utilise 
this interference effect to construct our entangled states and analyse in detail the impact of 
dissipation and detection efficiency on its stability. Our analysis of loss in this study goes 
well beyond that of previous works and allows us to realistically apply DIT to a quantum 
plasmonic scenario. 
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2. Physical system 

The physical system for the entanglement generation scheme consists of an array of spherical 
nanoparticles embedded in a dielectric material, with two distant QDs and four tapered metal 
nanowires, as depicted in a top-down view in figure [T] (a). All metallic regions have a 
frequency-dependent permittivity e m {co), and the background dielectric region has a static 
real and positive permittivity e d . The tapered metal nanowires are connected to adjacent metal 
nanoparticles and serve as input/output ports. The nanoparticle structure can be decomposed 
into two arms, each consisting of a linear array of nanoparticles, as considered in Ref. lfT2ll . 
which meet in the centre at a special arrangement of four nanoparticles that serves as a 
plasmonic beam splitter, as discussed in Ref. lfT3ll . Each arm of nanoparticles supports 
electron-charge density oscillations in the longitudinal and transverse directions with respect 
to the array orientation. For the orientation of the source nanotips shown in figure [T] (a), due 
to the direction in which the electron charge density oscillates in the nanowires, the nanowire 
field couples predominantly to the longitudinal oscillation in the nanoparticle. For coupling 
to the transverse polarisation, we rotate each source nanotip clockwise or counterclockwise 
by 90°. Each polarisation direction is maintained through the beam splitter whose corners are 
bent by 90° due to near field nature of the coupling with negligible radiation losses into the 
far field, provided that only nearest- neighbour interactions are considered (2Tl[22l[23]]. Thus, 
longitudinal (transverse) mode from the sources will become transverse (longitudinal) mode at 
the drains, and subsequently the drain nanotips should be properly oriented to collect charge- 
density oscillations in the corresponding direction. Alternatively, the out-of-plane transverse 
mode (perpendicular to the plane) can be used even if the corners at the nanoparticle beam 
splitter or the junction between the beam splitter and the two arms have arbitrary bending 
angles. The coupling strength between the nanoparticles depends on their separation distance 
and the direction of the electron-charge density oscillations HHEH. Each QD is coupled to 
the first metal nanoparticle of its respective array (as shown in figure [T] (a)), and assumed to 
have three internal energy states: a ground state, a long-lived metastable state, and an excited 
state, which we refer to as \m) and \e), respectively, as shown in figure [T](b). The states 
\g) and \m) represent the qubit states of the QD. The transition between \g) and \e) is assumed 
to be optically coupled to the adjacent nanoparticle, while the optical transitions connected 
to \m) are not. The desired level structure can be realised in a variety of solid-state material 
systems G5ll . 

2.7. Basic operation of the entanglement generation scheme 

If we take a single arm, say the left-hand side of figure [T] (a), then depending on the state of the 
QD, both the transmission and reflection coefficients of the metal nanoparticle array exhibit 
different characteristics. When the QD is in the state |m), i.e. when it is decoupled from the 
adjacent nanoparticle, the characteristics of the transmission and reflection coefficients of the 
source nanowire reveal a similar behaviour to those of a single array of nanoparticles (see 
Ref. (El), where the plasmonic field is transmitted across the array when it is on-resonance 
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Figure 1. (a) Setup for robust- to-loss entanglement generation using an array of metallic 
nanoparticles. The nanoparticle structure can be decomposed into two arms, each consisting 
of a linear array of nanoparticles, and a special arrangement of nanoparticles in the centre 
that serves as a plasmonic beam splitter. Two tapered metal nanowire waveguides on the left 
and right-hand side (source 1 and 2) focus light to the ends of their tips and excite localised 
surface plasmons on the adjacent nanoparticles. These excitations propagate across the arms 
and are then mixed at the nanoparticle beam splitter, after which they exit via the other two 
tapered metal nanowire waveguides (drain 1 and 2). By initialising the quantum dots (QDs) 
to be in a superposition of ground and excited states, and injecting coherent states of light 
into both the source nanowires, the system can be tailored so that a detection event at drain 1 
signals the generation of entanglement between the QDs. For the detection, an off-chip photon 
detector is linked to the output signal of drain 1, or alternatively this can be replaced with an 
on-chip detection unit l26l to directly measure the excitation at the (n + l)-th nanoparticle (as 
shown in panel (c)). A detection event at drain 2 (or alternatively on-chip detection linking to 
the {n + 2)-th nanoparticle) can also lead to entanglement generation between the QDs. All 
metallic regions have a frequency-dependent permittivity e m (oj), and background dielectric 
regions have static real and positive permittivity 6^, as defined in the text, (b) Energy level 
diagram for the QDs. Each QD (i = 1, 2) has three states \g t ), \mi) and \e t ) and has a \g)-\e) 
(\m)-\e)) transition with frequency (m) , which is detuned by (m) from the natural frequency 
of the adjacent nanoparticle, ojq. 



with one of the eigenmodes. On the other hand, when the QD is in the state i.e. when 
it is coupled to the adjacent nanoparticle, the propagating plasmons cannot pass through the 
adjacent nanoparticle and are completely reflected back into the source nanowire due to a 
phenomenon called dipole induced reflection (DIR), which is an alternative form of dipole 
induced transparency (DIT) EDI . Here, DIR occurs when the operating frequency oj of the 
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plasmons is resonant with the \g)-\e) transition frequency (as shown in | Appendix A[ l). The 
entanglement generation scheme then works as follows: Consider that light in the form of 
plasmon excitations are injected into the system via both of the source nanowires at the same 
time and that the QDs are prepared in a superposition, + N)i) ® ~^(lg)2 + 

If a plasmon is detected coming out from drain 1 for given input states at the sources, for 
example coherent states, then the plasmon came from the cases where either one of QDs is in 
|m), or both are in |m). A destructive interference of the case \m)\ ® |m) 2 can be induced by 
controlling the phases of the initial coherent states, so that only states \m)\ (8>|g>2 and \g)\ ®|m) 2 
lead to the detection event and consequently a superposed state of \g)\ ® |m) 2 and \m)\ <8> \g) 2 
is generated after the detection, i.e. \if/~) = -±=(\m)i ® \g) 2 - \g)\ ® |m) 2 ), which is an entangled 
state of the QDs. Likewise, the detection of a plasmon coming from drain 2 can be also 
used for generating the entangled state 1^"). This requires different phases of the initial states 
compared to the case of the detection at drain 1, and we focus on the case of the detection 
at drain 1. This is the basic concept for generating entanglement between two QDs using 
the plasmonic modes of a metal nanoparticle array in the ideal case when there is no loss. 
Here, we have made use of the characteristics of the transmission and reflection coefficients 
depending on the internal state of the QDs and an appropriate detection event. 



2.2. System Hamiltonian 

Before describing the entanglement generation scheme in more detail and showing its 
robustness to loss from the nanoparticles, we first briefly introduce the model to describe 
the system in figure [T] (a). We begin with the Hamiltonian for the nanoparticle system which 
is given by 

2n+2 

Hn P = J] ho)ia]ai + J] hgij (a]aj + a\a^j , (1) 

i=l [ij] 

where is the natural frequency of the field oscillation at the i-th nanoparticle, g t j is the 
coupling strength between the fields of the i-th and j-th nanoparticles, and [/, j] denotes a 
summation over adjacent neighbours j for a given nanoparticle i. The natural frequency ajt 
satisfies the Frohlich criterion Re[e m (<j)i)] = -2e d [HH|27]], which considers the nanoparticles 
to be small enough compared to the operating wavelength such that only dipole-active 
excitations are important [28J. The dielectric function of the metal e m (cxj) is given by a 
Drude-Sommerfeld model, which gives rise to a best fit to experimental data at frequencies 
corresponding to free space wavelengths /l > 350 (530) nm for silver (gold) ll29l . For 
simplicity we take all nanoparticles to have the same permittivity e m and radius /?, thus the 
local frequencies are set to be equal, cji = oj , for i = 1, • • • , 2n + 2. The operators a] (at) 
represent creation (annihilation) operators associated with a dipole-field excitation at the i-th 
nanoparticle, which obey bosonic commutation relations [a,-,at] = Sy. Here, a macroscopic 
quantisation of the fields is used, where the field modes are defined as localised solutions to 
Maxwell's equations satisfying the boundary conditions of the metal-dielectric interface [30J. 
In this case, the electron response in the metal is contained within the dielectric function of 
the metal, e m (aj) [EES I3T1 l32t 
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The interaction term in equation ([I]) involves two approximations for the center-to-center 
distance, d, between nanoparticles: a point-dipole approximation (37? < d), where multipolar 
interactions are negligible Il33lh and a near-field approximation, (d <c A), where X is the 
wavelength of the nanoparticle dipole field US [9), where the nearest- neighbour interaction is 
dominant via the Foster fields with a d~ 3 distance dependence [34]. These two approximations 
for d are covered simultaneously by a weak-coupling approximation, where the couplings 
between nanoparticles, g t j, are much less than the natural frequency, ojo, 0~H [T3ll . i.e., 
\gij «: (l>o for which we set max|g,-j| = O.Icjo throughout the work. 

The Hamiltonians for the QDs are then given by 

H 8 Q ^ = i»Qf (m) fr 8 J m) + hJ* (m) (fliof (m)t + a\&l (m) ) 

= \h&l + hJ\ « (a 2n+2 ^ ™ + al +2 a\ «) , 

where (m) is the transition frequency between |g) and |e) (|m) and |e» of QDi on the lefthand 
side,of (m) = \g)(e\ (\m){e\\&\^ = (|e><m|) and fr 8 ^ = k><e|-|g><g| (k)^|-|m)(m|). 
The vacuum Rabi frequency of QDi when coupled to the first nanoparticle, J\ (m) , is 
obtained within the rotating-wave approximation and the dipolar approximation, where the 
distance between the nanoparticle and QDi is assumed to be larger than the radius of the 
nanoparticle QUI . Finally, the dipole moment of the QD is aligned with the polarisation 
direction of the plasmonic excitation at the adjacent nanoparticle. Similar denotations and 
considerations as given above are used for QD 2 on the righthand side. 

2.3. Equations of motion for the system 

With the Hamiltonians introduced, the total system including the nanoparticle baths (for 
modelling loss) and the nanowires (sources and drains) can then be described in terms of 
Heisenberg equations of motion by using input-output formalism [jT2l [35l [36J] . The Heisenberg 
equations for the lefthand side of the system (from 1st to (n + l)-th nanoparticle and QDi) are 
given by 

— = -(ico + y + y)fli - igi,202 + VftTVi + V r i^in,i - iho\ (2) 

^ = -(io) + y )fl 7 - - igjj-iaj-t - ig jJ+1 a j+ i + ^jh^ (J = 2, 3, • • • , n - 1) (3) 
da^i § F^i . . / 

= -(iGJo + y)^ ~~ ign,n-l^n-l ~ ign,n+l^n+l ~ ign,n+3^n+3 + y^nhn,n (4) 

da n +\ . r^-fi /v . /v /v ~ / ~ 

— = -(/^o + y + — y )Sn+l ~~ ign,n+l^n ~ ign+l,n+2&n+2 + V&h^iiU + V^am- 1 ^in,n+ 1 (5) 

j - g (m) g (m) 

— ^- = -(in? (m) + A— )df (m) + Uiitid*™ + if (m) . (6) 

Similar Heisenberg equations can be written for the righthand side of the system (from (n + 2)- 
th to (2n + 2)-th nanoparticle and QD 2 ). The couplings between the nanowires and their 
adjacent nanoparticles are represented by g si , g S2 , g dl and g d2 . They are also assumed to 
obey the weak-coupling approximation, \g Sl \, \g S2 \, \g dl \ and \g d2 \ «: a) , as for the interparticle 
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coupling strengths g t j. This is imposed in our model for the nanowires, as the field profiles 
at the nanowire tips are similar in form to those of the nanoparticles, so that the model is 
a consistent description lfT2ll . The damping rate T for each metal nanoparticle depends on 
its size, and is given by Matthiessen's rule [TlTTl : Y = v F /A B + v F /R, where /l B is the bulk 
mean-free path of an electron, v F is the velocity at the Fermi surface, and 7? is the radius. 
By considering all the nanoparticles to have the same size and permitivity enables us to set 
= r for / = 1, • • • , 2n+2. The decay rate y g (m) of a QD from \e) to \g) (\e) to \m)) is given by 
the standard Wigner-Weisskopf spontaneous emission rate of the dipole, where it is assumed 
that the nanoparticle does not modify it. The operators ^ n ,i(2), ^in,i(2), and ^in,i represent 
the input fields of the source nanowires, the drain nanowires, and the /-th nanoparticle 's 
bath mode, respectively. The output fields of the nanowires, £ ut,i(2) and d ut,i(2)> and the 
nanoparticles' baths, £ out?i , are related to the input fields by a^n+i) = y==(Sin,i(2) + £ ut,i(2)X 

On+Kn+2) = -^=(^in,i(2) + ^out,i(2)), and a,- = -^(Vi +^out,i), for i = 1 2n + 2, respectively. 
Insertion and outcoupling losses in the nanowires are not included in the above Heisenberg 
equations, but will be included and discussed later. The damping rates of the QDs are assumed 
to be relatively small compared to all other rates, so that the noise operators, ff (m) (i = 1, 2), 
can be neglected in that they do not significantly affect the quantum coherence of the system. 

2.4. Scattering matrix for the system 

The entanglement generation scheme works in the weak excitation regime, where the dipole is 
assumed to be predominantly in the initial state \g) or \m). In this limit, o\ (m) can be replaced 
by its mean value, (erf (m) ) « -1, so that the coupled Heisenberg equations of motion give rise 
to four scattering matrices, one for each of the four internal states of the QDs (|mm), \gm) 9 
\mg) and \gg)\ with the general form 

Ci(^) = ^i Sl M^u U M + t SlS2 (aj)sl ut2 (aj) + t Sldl (aj)dl utl (aj) + t Sldl (aj)dl ut2 (aj) + f* uUsi (a>) 
sl n2 (a)) = t S2Sl (aj)sl utl (aj) + t S2S2 (o))sl ut2 (a)) + t S2dl (a))dl utJ (a)) + ^Mi^) + fL,^ *), 

(7) 

where the noise operators for the nanoparticles are j^ ut?Sl (^) = Zf* +2 ^sibi(^)^o Uti (^) and 
f^ uUS2 (w) = Ya n+1 4 2 bi(^)^outi^^ coefficient t D<> denotes the transition amplitude from an 
input mode "□" to an output mode "0". According to the state of the QDs |gm), |mg), and 
|mm)), the transition amplitude from "□" to "0" when the states of the two QDs are in Ix)®!}?), 
for ijg {m, g} is represented by the coefficient t^. For the sake of generality we have defined 
a dipole operator <j^ 2) for the dipole state |m), with the \m)-\e) transition frequency Q^ (2) and 
damping rate y^ (2) . However, in what follows we will treat the system when QDi (2 ) is in state 
\m) as if there is no QDi (2 ), since it is effectively decoupled from the nanoparticle in this case. 
For simplicity, we choose all the couplings between nanoparticles in the two arms to be equal, 
8j,k = gn P (at a fixed distance rf, array orientation, and nanoparticle size [fTTTO . and the nanowire 
to nanoparticle couplings, g Sl = g S2 = g dl = g S2 = gin-out- The weak-coupling approximations 
for these coupling strengths are equivalent to |g np | «: ojo and Igm-outl <^ <^o- We impose 
these by setting max|g np | = O.Iojq and max|g in _ out | = O.Icjq, which are achieved by varying 
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the distance between the nanoparticles, and between the nanowire tips and their adjacent 
nanoparticles, respectively, for given polarisations of electron-charge density oscillations lfT2ll . 
In this work we will mainly focus on the case of Aoj = co - co = 0, which is approximately 
the same as the case of using an optical pulse with a small enough bandwidth around co 
such that all coefficients in equation ([7]) are slowly varying. The transitions from the ground 
(metastable) states to the excited states for the QDs are detuned by 5\ (m) and 6 g 2 (m) from the 
resonant frequency coo of the nanoparticles, as shown in figure [T](b). 

From equation ([7]), the relation of the operators for the plasmons at resonance, co = co , 
can be written as a reduced matrix for each of the four internal states of the QDs to give 




^sidi(^o) 4,d 2 (^()) II -ouU v ~ u ' (g) 
*s 2 di(G>o) *s 2 d 2 (^o) /V ^Iut,2^ 



where k Sl(2) di(<^o)l 2 + ksi (2) d 2 (<^o)l 2 < 1 as the fields are reduced due to the losses from the 
metal nanoparticles and the QDs, as well as back reflections to the source nanowire modes 
■^out i(2)* The equality holds only if there is no dissipation and back reflections. A 50/50 beam 
splitter of lossy nanoparticles is achieved by setting g n , n+3 = g n+hn+2 = gh = (gin-out + r )/2 
and g n ,n+i = gn+2,n+3 = gv = V2g h , while taking care of the polarisation-dependence of the 
couplings between nanoparticles, such that t Sld2 (a>o) = ^ Sl di(^o) and f S2 di(^o) = ^s 2 d 2 (<^o) in 
equation ([8]) regardless of the internal states of the QDs. The coupling values for the 50/50 
beam splitter set the transition amplitudes t SlS2 (a>o) = t S2Sl (a>o) = in equation ([7]), which 
implies that any plasmonic field coming from the lefthand (righthand) side source nanowire is 
not influenced by the state of the QD on the righthand (lefthand) side, i.e. we have f* y = 

and C = *5- 

The role of the 50/50 beam splitter in our scheme is to erase the 'which- way' information 
of the two pathways associated with the QD internal states \gm) and \mg). Hence, the amount 
of energy that leaks into other output modes at the beam splitter is not important as long as we 
can perfectly erase the which- way information by having the same transition probabilities for 
the two pathways. Such a lossy beam splitter of nanoparticles is based on the lossless beam 
splitter proposed by Yurke and Kuang [13] . For the case Aco ± 0, the above choice of coupling 
values for g h and g v does not guarantee that the nanoparticle configuration performs as a 50/50 
beam splitter, i.e. it leads to a failure in the ability of the beamsplitter to erase the which- way 
information. This introduces decoherence into the generated entangled state. In figure [2} 
concentrating on the left arm and the beam splitter of nanoparticles, we plot the transmission 
and reflection spectral profiles for plasmons injected into the left source nanowire, \t 



\2 



sisil 



ks lS2 l 2 , ksidj 2 and |4 ld2 | 2 as Aaj/g np is varied, with gin-out/gnp = 0.5, r /g n p = 0.1, J\/g np = 0.3, 
and yi/gnp = 0.001 chosen as an example, when QDi is in \m) or \g). Here, we focus on 
the case of n = 2 and n = 3 in the left arm as the representative of even or odd number 
of nanoparticles because the dependence of the transmission and reflection amplitude on the 
in-out coupling strength gi n - ut is different for even and odd due to the resonances of the array 
(as shown in | Appendix A[ 2). For the case of QDi in state \m) (figure [2] (a) and (b)), on 
resonance (Aoj = 0) the left arm of nanoparticles follows the characteristics of a single array 
of nanoparticles and just transfers the plasmon field through the array to the beamsplitter 
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Figure 2. Transmission and reflection spectral profiles for plasmons that are injected into the 
left source nanowire. Left (right) column shows the case of n = 2 (3), and the upper (lower) 
row represents the case when QDi is in \m) (\g)). Regardless of the number of nanoparticles or 
the state of QDi, at resonance (Aoj = 0) the beam splitter splits the transmitted light field into 
two halves and blocks any field from reaching the right-hand side arm, so that | 2 = k Sl d 2 l 2 
and |4 lS2 | 2 = 0. In (a) and (b), when QDi is in |m), each coefficient exhibits the characteristics 
of a basic metal nanoparticle array, where n = 2 has |f SlSl | 2 while n = 3 has |f SlSl | 2 « (see 
main text for details). In (c) and (d), when QDi is in \g), DIR causes \t SlSl | 2 « 1 for both n-2 
and n = 3 when Aoj = 0. 



where it is split equally. The reason for the large backscattering (|f SlSl | 2 probability) in the case 
of n = 2 (figure [2] (a)) is due to the eigenmodes of the array (from the source to each drain) not 
being on resonance for Aoj = 0, whereas for n = 3 (figure [2] (b)) there is a resonant eigenmode 
at Aoj = for each drain lfT2l . For Aoj ± the various eigenmodes of the system come into 
resonance at particular frequencies of the input field and cause the oscillatory behaviour seen 
on either side of Aoj = 0. On the other hand, for the case of QDi in state \g) (figure [2] (c) 
and (d)), the first nanoparticle is coupled to QDi and DIR occurs (as shown in | Appendix A| ), 
regardless of the resonant eigenmodes in the nanoparticle system. Thus in both figure [2] (c) 
and (d) at Aoj = the backscattered probability |£ SlSl | 2 goes to unity. Away from Aoj = the 
transition probabilities follow the same behaviour as for the case of QDi in state |m>, as the 
QD becomes off resonant and so cannot induce any dipole-based interference effect. 



3. Entanglement generation scheme 



In this section, we first consider a weak input classical source to delineate the basic idea of 
how the scheme works and then describe a more general scenario. 
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3.1. Weak field injection 

To generate entanglement between the QDs, both are initialised to be in an equal superposition 
of qubit states \g) and |m), i.e. |+) = ~^(\g) + N». This can be achieved by first driving the 
QDs into the lowest-energy state and then rotating them by either a direct n/2 transition, or a 
Raman transition E71 . depending on the specifics of the QD internal energy level structure. 
After initialising the QDs, a coherent field \a) with frequency ojo is injected into the left source, 
^ nl , and simultaneously another coherent field \/3) with same polarisation and frequency as 
\a) is injected into the right source, s^ T The initial state of the whole system is then 

m = \(\g)i + \mh)(\gh + Im^lar^^lO^IO)*. 

The input plasmon fields interact with the QDs only if they are in the state \g). In the case that 
one (or both) of the QDs is in the state |m), a field is transferred along one (or both) of the 
arms, then mixed at the beam splitter and goes out through one of the drain nanowires, d^ ut 1(2) . 
Entanglement can be generated by the detection of an excitation in drain 1, which could have 
come from either the left or right arm, when the QDs were in \mg) and \gm) respectively (the 
case when both arms supply an excitation, with the QDs in |mm), is suppressed as described 
below). Thus, the nanowire mode ^ ut x is connected to a detector in order to postselect the 
successful case when the which-way information has been erased and the QDs have been 
entangled as a consequence. Alternatively, an on-chip detector could be directly coupled 
to the (n + l)-th nanoparticle so that the drain nanowire is not needed. Note that the total 
plasmonic system only requires a classical source of light in order to generate entanglement. 

To see how the scheme works, it is useful to consider sufficiently weak coherent states 
(a «: l,yS «: 1) when there is no loss present. We may expand |^) to first order in the 
plasmonic excitation number and approximately drop all other output modes. Then, as the 
scheme is based on the detection of an excitation we can drop the vacuum terms so that 
|Y/> * l(\gh + |ro>i)(|g>2 + \m) 2 ){a\\) Sl +/?|1>, 2 ). Perfect DIR leads to t 88 ^ = t 8 ^ = and 
*sfdi = CI = 0> an d using equation ^ for each of the internal states of the QDs, together with 
a postselection of the state when an excitation is present in drain 1, we have the unnormalised 
state of the QDs as 

I^qds = \{{a^ +PQ i )\mm) + /?Qgm> + 04 (9) 

Here, the beam splitter enables t™™^ = t™ 8 d and = ^ , i.e. the plasmon injected in 
each source is not influenced by the QD in the opposite arm. If the phase of the amplitude 
yS is chosen such that af^ + Pf^l = 0, we have a normalised ideal entangled state as 
I^)qds = ~fi(\ m 8) ~ \§ m )) - l*A~) w ^h probability |a^J 2 . This is the basic idea of the 
scheme for the limiting case of weak a and /?. 

3.2. Arbitrary field injection and loss 

We now turn to the more general scenario, starting with the initial state |^-) without any 
restrictions on a and J3. In this case the detection of excitations at drain 1 is modelled by 
the projection operator P dl = J^=i l n )di(w| = I - |0> dl <0| that projects the state of the system 
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onto a subspace containing at least one excitation in mode dl uiV This projection models the 
measurement performed by an ideal non-photon number resolving detector, which registers 
a detection event as long as there is at least one excitation at the detector. Later we will 
also investigate the effects of detection inefficiency in the scheme. As a result of lifting the 
restriction on the input field amplitudes and the accuracy of the detection, the state of the QDs 
becomes mixed and can be written as 

11 (all fields) 

Tr(QDs)[Tr( a n fie i ds ) 

The final state of the output fields, is found by transforming |^) according to the four 
sets of scattering matrices from equation ([7]) depending on the state of the QDs. It is written 

as 

P*7> = \(\ggM 88 hzUs + Nm)ir m )fields + lgm>|^>fields + Ng> ^fields), (11) 

where the state of the fields at the output modes is represented by a product of coherent states 

Infields = |^) Sl l^>s 2 l//f>d 1 l//f>d 2 ®Jf \XJ\, (12) 

with amplitudes = at^ s . + ]3t^ s . (j = 1,2) at the source nanowires, fij = at^ & + fit^ & 
(j = 1, 2) at the drain nanowires, and^J = at^ h _ +/fr^ b . = 1, • • • , 2n + 2) for the nanoparticle 
bath modes. Each coherent field consists of an interference of two pathways, one from source 
1 and the other from source 2 with respective transition amplitudes. Throughout this work, we 
impose the matching condition <rf™ + f}t™£ = in order to exclude the possibility of the case 
of \mm) leading to a detection, i.e. the amplitude for the coherent state in drain 1 is jji™ m = 0. 
As mentioned previously, this can be achieved by properly adjusting the amplitude ft of the 
coherent field injected into source 2. In order to measure how well the QDs are entangled, 
we employ the fidelity as a figure of merit, defined as the overlap between the desired final 
state, \if/~) 9 and the actual final state of the system, Pqd s , and given by F = (^"Ipqds l^ - )- 
The fidelity allows us to quantify how close the final state is to the desired state, and at the 
same time provides a lower bound on the concurrence, C, a typical entanglement measure 
for two qubits, i.e. C > max(0, 2F - 1). Such a lower bound implies that entanglement 
is always found when the fidelity is greater than 1/2 [[38l [391 . The fidelity is given as 
F = <<A"IpqdsI*A"> = ^(P22 + P33 - P23 - P32), where p pq is the entry in the p-th row and 
the q-th column of the matrix for p QDs that is spanned by the basis {\gg) 9 \gm) 9 \mg) 9 \mm)}. 
The entries are given by 

P22 = ^(l - ^riO)d 1 <OK m )d 1 ) 5 P33 = ^(l - <Pri0> dl <0K%) (13) 

P23 = ^(KX'V - </4 m io> dl <oK*> dl ) x 

The efficiency (or success probability for generating an entangled state) is the probability to 
detect photons at drain 1 and is given by the denominator of equation ( fT0| ) as 

7] = Tr ( Q Ds )[Tr (a ii fie i ds) [Pd 1 l^ / /)(^ / /|PdJ] 

= l - U\(0\txf )| 2 + kok™>I 2 + KOIaO 2 + l<0K*>l 2 ]. (14) 
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Figure 3. Fidelity of the state of the QDs, pqd s , with respect to the maximally entangled 
state at resonance Aoj = for the coherent fields injected into the sources. Here, the 
magnitude of the field of the injected states, a (with j3 = ia), is varied for n = 2 and n = 3. 
As an example, we have used the system parameters gin-out /gnp = 0.5, Jf/gnp = J\lgnp = 0.3, 
yf/gnp = yf/gnp = 0.001 and To/gnp = 0.1. When a is increased, the fidelities drop and 
asymptotically approach 0.5, whereas the efficiencies go up and asymptotically approach 
0.5. The gradients of the fidelities and efficiencies can be changed depending on the system 
parameters we have set. These affect the performance of the scheme, as shown in figure [4] 
where we vary g in -out/gnp and r /g n p. 



Thus, by inspecting equation ( [13] ), a high fidelity is obtained in the limiting scenario where 
we have /u 8 ™ = -/i™ 8 with all the other field modes for the case when the QDs are in \gm) 
approximately equal to the field modes for the case when the QDs are in \mg), so that they are 



effectively factored out in equation ( pTj ). In this case, the state of the total system becomes 
|*F/) = (/iflgm) + fi™ 8 \mg)) ® |all fields). However, this limiting scenario is hard to achieve 
for arbitrary fields injected into the system as higher-order photon number contributions act 
as a decoherence mechanism on the final state. The impact of higher-order photon number 
contributions on the value of the fidelity and efficiency are shown figure |3j Here, a trade- 
off between fidelity and efficiency can be seen as a is increased (note: ft = ia by the 
matching condition af™l + fif^ x = 0) BOll . The fidelity decays quickly with increasing a 
due to the presence of higher-order excitations in the plasmonic system. It asymptotically 
approaches 0.5 in the limit a » 1, indicating that the higher-order excitation contributions 
have completely decohered the state of the QDs. 

4. Results 

Having described the entanglement generation scheme and general trends for the fidelity 
and efficiency we now go into the details of the robustness of the fidelity to loss from the 
metal nanoparticles, the length of the arms, and inhomogneous broadenings of the QDs. 
We also examine the efficiency of generating entanglement in these circumstances. Our 
analysis of entanglement generation is classified into two groups, each according to the 
number of nanoparticles in the arms: even or odd. This is because the dependence of the 
transmission amplitude on the in-out coupling strength gi n - ut is different for even and odd 
due to the resonances of the array (as shown in | Appendix A[ 2). We have also chosen to 



plot all parameters in units of the nanoparticle coupling g np so that the plots are independent 
of g np , as long as g np «: ojq is satisfied (weak-coupling approximation). Note also that 
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Figure 4. Fidelity of generated QD state, pqd s , with respect to the maximally entangled state 
when the input fields are on resonance (Aoj = 0) as the in-out coupling strength gin-out/gnp 
and the amount of metal loss r /g n p are varied, for n = 2 and n = 3. Left (right) hand column 
shows the case of n = 2 (3), and the upper (lower) row shows the fidelity (efficiency). 



max(g in _ out /g™ ax )=l must be imposed, where g™ ax = O.Icjo, otherwise we move away from 
the weak-coupling regime for the sources and drains. In other words, the rescaled couplings 
gin-out/gnp can in principle go higher than 1, but the value for g np must be lower than O.lcoo 
to compensate, so that we are still in the weak-coupling regime. All calculations are done for 
resonant input fields, Aoj = co-co^ = 0, except in subsection |4 . 5 1 where we investigate the pulse 
profiles of the inputs. We also set the QD-nanoparticle couplings as /f/g np = J\ /gnp = 0.3, 
and the QD decay rates y\ /g np = y g 2 /gnp = 0.001. Other parameters of the system are varied 
in our analysis. 



4.1. Robust fidelity against metal loss 

We first examine the fidelity and efficiency of our scheme with respect to loss in the metal 
nanoparticles. For this purpose, we choose the cases oin-2 and n - 3 as representatives of 
even and odd numbers of nanoparticles in each arm. Figure [4] shows the fidelity and efficiency 
as the nanoparticle loss r /g np increases and the in-out coupling strength gin-out/gnp is varied 
for a = 0.5 (with {3 = ia). Here, we have assumed that the \g)-\e) transition of each QD is 
resonant with its adjacent nanoparticle such that #[ (2) = 0, while 5^ (2) is chosen for the QD 
\m)-\e) transition to be completely decoupled from the dynamics. 

In figure [4] (a) and (b), we show that the fidelity varies only slightly with increasing 
To/grip, implying a remarkable robustness against metal loss. Such a robustness in fidelity can 
be understood from the observation that the reflection coefficient of each arm predominantly 
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determines the fidelity for the case of resonant QDs (6\ = 6 8 2 = 0), as we shall explain in 
detail below. This indicates that the fidelity mostly depends on the behaviour of the reflection 
coefficient of each arm with metal loss T /g np and gin-out/gn P ; the reflection coefficient varies 
only slightly with T /g np since the energy never enters the array, and has a monotonic trend 
with gin-out/gnp (see | Appendix A| ), resulting in robustness of the fidelity against metal loss and 
a general decrease in the fidelity with respect to gin-out /gnp> respectively. Contrary to the case 
of the fidelity, we show that the efficiency of generating entanglement is decreased with r /g np 
in figure [4] (c) and (d). This is because the efficiency is determined by how much energy is 
transferred from the sources to the drains, indicating that the efficiency is mostly related to the 
transmission coefficient of each arm. In contrast to the reflection, the transmission coefficient 
is more sensitive to metal loss, and has different dependence on gm-om/gnp (as shown in more 
detail in |AppendixA| ) . 

To understand the above observations in detail, we analyse the final state given in 
equation ( 11 ). First, let us consider the fidelity. For the case of resonant QDs (6 8 x = 5 8 2 = 0), 
we have the coherent amplitudes for drain 1 as fi 8m - -fi™ g , a condition required to generate 
a perfect singlet state, \if/~), as seen in the discussion of the 'limiting scenario' in the previous 
section. We also have the coherent amplitudes for drain 2 as fi 8 ™ = fx™ 8 , so that the 
corresponding term in equation (fTTl) can be factored out from On the other hand, the 
coherent amplitudes for the other output modes, corresponding to the sources, £ gm (= -it;™ 8 ) 
and £™ g (= -i^l™), are not equal to each other and cannot be factored out, which leads 
to dephasing of the QD state. In addition, all of the bath modes (modelling loss at the 
nanoparticles) are also involved in dephasing the QD state in a similar way. These bath modes 
dominate the value of the fidelity for large values of r /g n p as p 2 3 (= p* 32 ) — > 0. Nevertheless, 
the dephasing effects of the bath modes and their impact on the fidelity are relatively weak 
compared to the source modes until large losses are incurred in the system, i.e. for r /g np > 1. 
For these reasons, both the coherent amplitudes for the source output modes, % 8 ™ and £ mg , are 
predominantly responsible for the behaviour of the fidelity before T /g np becomes very large. 
Thus, we focus on investigating the behaviour of the amplitudes £ gm and £ mg . Due to the beam 
splitter, the transition amplitudes *f 2 ™ = and t™ 8 Sl = 0, so that = at 8 ^ and £™ g = art™, 
where t 8 ™ Sl is the reflection coefficient when DIR occurs for QDi, while t™ 8 Sl is related to the 
reflection coefficient of a single arm in the absence of QDi. Thus, the more similar t™ 8 Sl and 
t%™ Sl are to each other, the more similar % 8 ™ and £ mg are to each other, and the output states in 
the corresponding modes can be factored out from the state |*F/) in equation (fTTl). This leads 
to a high fidelity. Otherwise, the more t™ 8 Sl and t%™ Sl are dissimilar to each other, the more they 
cause dephasing. 

To compare £ 8m and % mg \ we choose as an example a = 0.5 in figure [5j and show the 
dependence of t 8 ™ Sl and t™ 8 Sl on the input/output coupling gin-out /§n P for a set of values of 
metal loss. In figure [5] (a) and (b), the coefficient (dashed line), which is related to the 
occurrence of DIR, does not change appreciably with respect to the increase in loss. This is 
because the first nanoparticle is not excited due to quantum interference when DIR occurs 
(see |Appendix A[ 2), similar to a waveguide coupled to a cavity with a dipole HDll . On the 
other hand, the coefficient related to the absence of a QD, t™ 8 Sl , exhibits the characteristics of a 
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Figure 5. Origin of the robustness of the fidelity. In (a) and (b) the transition amplitudes t™ 8 Sl 
and ff^ are shown. In (c) and (d) the density matrix entry p 2 3 is shown. In all plots gin-out /gnp 
is varied for set amounts of metal loss: To/g np = 0.01, 0.1, 0.5, and 1. 



nanoparticle array, where the dependence on gin-out /g np is different for even and odd numbers 
of nanoparticles, as mentioned in the previous section. For n = 2, the coefficient t™ 8 Sl diverges 
from 4™ Sl as gin-out/gnp increases and reveals a slight variation in r /g n p, which explains the 
falling fidelity as gin-out/gnp increases and robustness of the fidelity with increasing r /gnp, as 
shown in figure [4] (a). For n = 3, the coefficient t™ 8 Sl also diverges from t 8 ™ Sl as gin-out/gnp 
increases, which again explains the fidelity decay as gin-out/gnp increases in figure [4] (b). 
However, contrary to the case of n = 2, the coefficient t™ 8 Sl moves closer to 4m as fVgnp 
increases when gin-out/gnp is small, and shows a bigger variation compared to the n = 2 
case when gi n -out/gn P is increased, as shown in figure [5] (b). This analysis explains the slight 
increase in fidelity with r /g np for small gin-out /gnp, and the robust fidelity with r /g n p for large 
gin-om/gnp that is shown in figure [4] (b). While £ 8m and £™ 8 are approximately responsible for 
the fidelity when 5\ = 5 8 2 = 0, the entry p 2 3 in the QD density matrix is the only entry that 
directly determines the fidelity among the different entries of Pqd s > since p 2 i = P33 = 1/2 
when 5\ = 8\ - 0. In figure [5] (c) and (d), we plot p 2 3. By comparing the behaviour of p 2 3 in 
these plots with that of the fidelity in figure [4] (a) and (b), it can be seen that as the real part of 
P23 goes to -i the fidelity increases. Thus, the behaviour of p 2 3 follows the trend of *™ with 
increasing r /g n p and gin-out /g np , except for the intermediate regime of gin-out /gn P - Here, there 
is an interplay of other field modes in the system, which become significant and make a small 
difference in the trend of p 2 3 compared to that of the source coherent amplitudes £ 8m and %™ 8 '. 

Secondly, let us consider the efficiency. In contrast to the fidelity, the efficiency depends 
only on the amplitudes jj 88 , \I[ m , jj 8 ™, and fi™ 8 , as seen in equation (14), which are mostly 
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related to the transmission amplitudes of the two arms. In |Appendix A[ 2, the square of 
the transmission amplitudes are presented in figure |A2[ where for n = 2 the transmission 
is found to have its highest value at gin-out /gnp = 2, whereas for n = 3 the transmission 
quickly increases with increasing gm-out/gnp- Such characteristics of the nanoparticle array 
are reflected in the efficiencies for both n = 2 and n = 3, as shown in figure [4] (c) and (d). 
Higher values of the transmission as gm-out/gnp is varied enable the probability of the injected 
plasmons to exit via drain nanowires to become more likely, so that the efficiency to detect 
excitations at drain 1 increases (compare figure [4] (c) and (d) with figure |A2| (a) and (b)). This 
shows that the efficiency does not have simple trade-off with the fidelity in figure [4} 

In this subsection, we chose a = 0.5 as an example. However, if we choose a larger 
value of a then higher-order excitations will reduce the fidelity and increase the efficiency 
on average in figure |4j as already shown in the previous subsection. Note that this is not a 
practical problem as the amplitude a is easily controllable, so that it can be decreased as low 
as required (or increased) for a reasonable fidelity and efficiency. Thus, using coherent states 
as initial states is a key merit of the scheme, which allows us to use a metal nanoparticle array 
supporting lossy plasmonic modes. 



4.2. Robust fidelity against the length of a nanoparticle array 

In terms of the total amount of loss present in the entire system, the effect of increasing the 
loss of each nanoparticle, for a fixed number of nanoparticles, can be regarded as equivalent 
to increasing the number of nanoparticles in each arm for a fixed amount of loss. In figure [6} 
for a fixed amount of metal loss, T /g np = 0.1, and input amplitude a = 0.5, the fidelity 
and efficiency are shown as gi n - ut /g np is varied with increasing number of nanoparticles. As 
expected, the behaviour of the fidelity and the efficiency as the length of nanoparticle array is 
increased is very similar to the behaviour seen in figure |4} It is quite remarkable that one can 
achieve a robust high fidelity (and thus entanglement) over a reasonably long array of metal 
nanoparticles, even though it has so far been believed that loss limits the length of an array of 
nanoparticles for quantum information processing purposes. 



4.3. Robust fidelity against detunings of the QDs 

A major challenge when using solid-state emitters is inhomogeneous broadening, typically 
caused by emitter size variation and strain fields in the host material. This means that 
two emitters usually have non-identical emission wavelengths and thus the two QDs in our 
scheme cannot be easily assumed to have the same resonant transition frequencies in a realistic 
scenario. Here, as one of the merits of the scheme, we show the robustness of fidelity against 
detunings of the QDs. Such robustness has been pointed out already in Ref. Il40ll . but here we 
provide a more detailed analysis for the case of a plasmonic nanoparticle array system. In this 
section, we set 5\ = 5 + A5 and 6 8 2 = 5 ~ A£, and assume that y\ - y 2 and J\ = J 2 . In figure [7] 
the fidelities are shown as the average of the detunings, 6 = (#J + 6 8 2 )/2, and the difference 
in the detunings, A6 = (6 g x - S g 2 )/2, are varied for a = 0.5. Here, we have taken the matching 
condition used in the previous section, or^J + /^di = ®' N° te ' however, that this matching 
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Figure 6. Fidelity of generated QD state, pqd s , with respect to the maximally entangled 
state \iff~) when the input fields are on resonance (Aa) = 0) as the in-out coupling strength 
gin-out /gnp is varied and the number of nanoparticles n is increased for a set amount of metal 
loss, To/gnp = 0.1. Left (right) column presents the case of even (odd) n, and the upper (lower) 
row presents the fidelity and efficiency respectively. 



condition does not guarantee the optimal fidelity. For example, in the limiting case of weak 
coherent states, dropping all other field modes, the final state given in equation (fTTl) can 
be rewritten as |*F/)qd s ~ \(l Ji \ g \gg) + fi™ m \mm) + \i g ™\gm) + fi™ 8 \mg)), which is similar to 
equation ([9]). For | v F/)qd s to become |t/r~>, we need to make sure that there can be no detection 
at drain 1 when both the QDs are in |m), or both are in \g). The matching condition implies 
there is no detection for the case of both QDs in state |m), i.e. jji™ m = 0, whereas for the case of 
both QDs in state \g) the probability of a detection at drain 1 can be removed by perfect DIR, 
i.e. = t^ d = — > fif = 0. If imperfect DIR occurs due to the detunings of the QDs, we 
can no longer guarantee that there will be no detection at drain 1 for the QDs both in state \g) 9 
i.e. fif •$> 0. For this reason, the matching condition does not guarantee the optimal fidelity, 
but it is still sufficiently useful to observe the robustness of the fidelity against the detunings 
of the QDs. An optimization of the amplitude /? to reach a higher fidelity than simply using 
the matching condition might lead to more robustness of the fidelity against the detunings of 
the QDs, but this involves numerical calculations that make it complicated to understand the 
origin of the robustness. Thus we do not consider it in this work. 

In order to understand the behaviour of the fidelities as the detunings vary, as shown in 
figure[7](a) and (b), it is helpful to analyse two extremal cases. The first is when the respective 
detunings of the two QDs have different signs while their average is zero, i.e. So/g np = and 
A5/g np + 0, labelled as region A (see figure [7] (c) for level diagram), corresponding to a cut 
along the y-axis. The second is when the respective detunings of the two QDs are the same 
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Figure 7. Fidelity of generated QD state, Pqd s , with respect to the maximally entangled state 
\\ff~) when the QDs are detuned. In (a) and (b) the input fields are on resonance (Aco = 0) 
as the average detuning Sq /g np and the difference in detuning AS /g np of the QDs is varied for 
n -2 and n - 3, respectively. Two extreme regions are labelled: A is when So/g np = and 
A£/gnp ^ 0, and B is when So/g np ^ and AS/g np = 0. The corresponding energy level 
diagrams are shown in (c). In (d) and (e) the fidelities in the regions A and B are shown, 
respectively, for both n = 2 and n = 3. In (f) and (g) the reason why the fidelity is more 
robust against a change in A£ /g np the more detuned £o /g np is from coo can be seen from the 
behaviour of the transmission and reflection spectrums for two separate values of £o /g np (solid 
and dashed) and a fixed value of A£ /g np 



as each other, i.e. 6 /g np ^ and A6/g np = 0, labelled as region B (see figure [7] (c) for level 
diagram), corresponding to a cut along the x-axis. 

In region A, imperfect DIRs at both ends of the array occurs due to the QD detunings. 
The asymmetric detuning leads to fif ^ 0, so that |^/)qd s ~ \(Mi g \gg) + ft 8 i m \g m ) + ft^lmg)) 
and there exists a non-zero probability of detection at drain 1 when both the QDs are in \g). 
Thus, the state \gg) gives a detection, causing a loss of fidelity. Furthermore, in region A we 
have fi 8 ™ = -(ji^ g y. This is not a problem if the respective imaginary terms are approximately 
zero. However, if this is not the case then it leads to dephasing in the final QD state. Here, the 
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phase difference between fi 8m and fi™ 8 affects the relative phase cp of the generated entangled 
QD state as -h^(\gm) + e l(f) \mg)). From this point-of-view, \a 8 ™ - -fi™ 8 is highly desirable for 
achieving a high-fidelity, as pointed out in the limiting scenario in the previous subsection. 
Thus, increasing |A5/g np | causes a loss of fidelity since the respective imaginary terms of fi 8m 
and fi™ 8 are increased, as can be seen in figure [7] (d). On the other hand, in region B the equal 
detunings of the QDs enables jj 88 = to be obtained via destructive interference with the help 
of the matching condition, i.e. |^/)qd s ~ ^iM^lgm) + ^^S))- It also enables /i 8m = to 
be obtained, which provides a high fidelity, in contrast to region A. Here, increasing |5o/gn P l 
disturbs the DIRs at both ends of the array, so that the amplitudes related to the occurrence 
of DIR are no longer much different from the amplitudes that follow the characteristics of the 
system in the absence of the QDs. Thus, all the other field modes for \gm) are approximately 
equal to those for \mg) and can be factored out so that \¥f) — > (ji 8m \gm)+ii™ 8 \mg))<g)\all fields). 
In addition, the magnitudes of fi 8m and fi™ 8 decrease with increasing |5o/gn P l such that 
moves closer to the state |^/)qd s obtained in the limiting scenario of weak coherent states. As 
a result, the fidelities increase with increasing |5o/gn P l for both n = 2 and n = 3, as shown in 
figure[7](e). 

Interestingly, the robustness of the fidelity against an increase in \A6/g np \ observed around 
region B is as follows; the greater |<Wgn P l becomes, the more robust the fidelity is against 
changes in \A6/g np \. This robustness can be understood by looking at how sensitive the various 
amplitudes are as they change from their values that lead to the best fidelity in region B when 
we vary A6. Consider two examples, So/g np = 0.02 and <Wgn P = 0.12, for which the change 
in the amplitudes as Aaj/g np varies is shown in figure [7] (f) and (g), for both n = 2 and n = 3. 
Here, it can be seen that by varying <Wgn P (solid and dashed lines represent two different 
values) the resonant peaks and dips in the transmission and reflection spectrums are shifted. 
In addition, for a fixed value of A6/g np for the two chosen detunings of 6 /g np , a variation in 
the spectrums is induced. As we are interested in the difference in the detunings, Aoj, one can 
see the values of the specturms at Aaj/g np - are changed less sensitively by a slight variation 
of A6/g np for 6o/g np = 0.12 than for 6o/g np = 0.02. Such characteristics of transmission and 
reflection spectrums result in the robustness of the fidelity against the QD detunings A6/g np . 

4.4. Nanowire loss and detection efficiency 

As we are interested in the mapping of the input fields at the source tips to the output 
fields at the drain tips in the scheme, we have assumed so far that the source and drain 
excitations experience no loss when propagating into and out of the tip regions. However, 
in a realistic system, the losses in the nanowires would affect the overall performance. Here, 
we consider two kinds of losses in the nanowires: (i) insertion loss, which occurs when the 
light is transferred from a far-field source into the source nanowires and the resulting surface 
plasmons propagate along the nanowires up to the 1-st or (2n + 2)-th nanoparticle, and (ii) 
outcoupling loss, which occurs when surface plasmons propagate along the drain 1 nanowire 
and are extracted into the far-field to be detected (or equivalently the detection efficiency of an 
on-chip detector directly connected to (n + l)-th nanoparticle). For insertion loss, this is not 
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(a) Fidelity for n=even (b) Fidelity for n=odd 




(c) Efficiency for n=even (d) Efficiency for n=odd 




Figure 8. The fidelities and efficiencies for both even and odd numbers of nanoparticles are 
shown as the detection efficiency (or nanowire loss), k, varies for a = 0.5 and gin-out /gnp = 0.5. 
The left (right) column shows the case of n even (odd), and the upper (lower) rows show the 
fidelity (efficiency). 



a problem for the scheme as we can simply compensate it by increasing the intensity of the 
input coherent states as much as we need for the plasmon excitations at the 1-st and (2n + 2)- 
th nanoparticles. In this sense, one of the advantages of the scheme is the use of coherent 
states as initial 'seed' states for the generation of entanglement. These enable the scheme to 
be stable no matter how much loss the initial coherent fields experience on input to the array 
system. On the other hand the second type of loss, the outcoupling loss, more significantly 
affects the fidelity and the efficiency. The outcoupling loss can be included as a 'lumped 
detection efficiency', k, by replacing P dl with P dl , where P dl = I - Z^ (l - K) n \n) dl (nl which 
corresponds to the case of ^-photons being detected by the non-photon number resolving 
detectors with a probability 1 - (1 - k) u Il39ll . Such a treatment of outcoupling loss is valid 
due to the fact that for an output coherent field at drain 1, say the expectation value of the 
measurement of an ideal detection with outcoupling loss amplitude is equivalent to that 
of a lumped detection without outcoupling loss, i.e. (yu|P dl |yu) = 1 - e~ K ^ = ( V^/|P dl | V*"A*>- 
This compatibility between nanowire loss and detection efficiency allows us to treat them in 
the same manner. To show how much the output coupling loss (or detection efficiency) affects 
the entanglement generation scheme, in figure [8] we plot the fidelity and the efficiency as the 
detection efficiency, k, is varied and the number of nanoparticles, n, increases for a = 0.5 
and gin-out/gnp = 0.5. As the detection efficiency k increases, the fidelity and the efficiency 
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become generally higher. The gradients of the fidelity and the efficiency as k varies depend on 
the magnitude of the input amplitude, a. For very small a (or large a) the gradient as k varies 
is small, whereas for a having an intermediate value, the gradient is big. 

4.5. Pulse width of input coherent fields 

We now investigate the effects of the pulse width of the coherent field injected into the 
source nanowires. We start by focusing on the coherent field that comes from the source 
nanowire on the left-hand side, whose multi-mode coherent state can be described by 
\{a}) Sl = expCsJ^ - Si n ,J|0) Sl , where the wavepacket operators are s\ na = JT daja(aj)sl n (a)), 
with J ^ d(^|a(a;)| 2 = (h a ) [41 j. For concreteness, we consider a Gaussian wave packet with 
spectral amplitude profile a(aj) = (2^cr 2 )" 1/4 ^" (wo_w)2/4cr % where ojq is the central frequency 
and cr a = 6co a /(2 V21n2) is the standard deviation corresponding to a full width at half 
maximum bandwidth 5co a for the spectral intensity profile |ar(<^)| 2 . 

Such a coherent field pulse should be carefully employed with appropriate constraints for 
the analysis of this work to be valid. First, for the monochromatic approximation (Aoj = 0) 
to be valid, 5co a should be narrow enough such that as the entanglement generation scheme is 
carried out, each amplitude of the system hardly changes (or varies slowly). This is equivalent 
to non-dispersive transfer of the plasmons in the metal nanoparticle array lfT2ll . Second, the 
higher-order photon number contributions that induce decoherence impose a constraint of 
(h a ) «: 1 for a high-fidelity entangled state to be generated between the QDs. Third, the 
analysis of the protocol using DIR is valid in the weak excitation regime where the QDs are 
unsaturated, i.e. (dj.(t)) « -1 which is equivalent to « 1. If the amplitude of the 

input coherent field is large enough such that the g-e transitions of the QDs are saturated, 
then they will lead to an optical nonlinearity and linewidth broadening ll42l . and equation ([6]) 
can no longer be treated in the linear regime. To investigate the implications of this weak 
excitation limit, we use the Heisenberg equations of motion for the first nanoparticle and QDi, 
and ignore the field operator for the second nanoparticle since it never has the chance to be 
excited when DIR occurs. Thus, equation ([2]) and ([6]) can be rewritten in the weak excitation 



regime as 



da i 

dt 
daj 




(15) 



_i = -(^ + ^) ( xf-// 1 a 1 . (16) 

Eliminating a\ from the above equations when (%^t + « we have at resonance 



(Aco « 0), 



- i5\(^Y^ + y) - J 2 \a\ = -Ui V<?in-outSin,i - iJi VfoVi. (17) 



The above equation can be multiplied by its conjugate to give 
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where {b[ n x b. m x ) « 0, {s\ n { b in { ) « 0, and 1 ^ in x ) ~ are assumed, and (& 8 ^&\) represents 
the probability of QDi being in the excited state. Here, {s\ n x s in x ) is identified as the total flux 
of photons in the input field \a) Sl of frequency aj . In the weak excitation limit, «: 1, 

which finally puts a limit on the excitation number (h a ) of 

J\ ^(gm-out + r ) 2 \a\ 2 (h a ) 

» 

4Jfg in _ 0Ut 6r a 6r a 



+ ^/rr u/ » ^ ^ ^r-, d9) 



where 5r a is the pulse width of |{^}) Sl . Note that this third constraint of the weak field 
excitation approximation can also cover the first and second constraints of sufficiently narrow 
bandwidth 6aj a (or equivalently long pulse width 6r a ) and low excitation number (h a ) «: 1. 
Similar constraints on {hp) for \{/3}) S2 can be derived from the Heisenberg equations of motion 
for the (n + 2)-th nanoparticle and QD 2 . 



4.6. Initialisation of the QDs 

In the scheme, the QDs are initially prepared in an equal superposition of \g) and \m). 
However, here we show that an equal superposition is not a strict requirement for the QD 
initialisation. When the QDs are prepared in arbitrary states, the initial state of the system is 
given by 

l^'i) = (c g i\g)i + c ml \m)i)(c g2 \gh + c m2 |m> 2 )|or>^ |y0>^ 2 10>^ |0>^ 2 , (20) 

where \c g i\ 2 + \c m i\ 2 = 1 and \c g2 \ 2 + |c m2 | 2 = 1. In the ideal case of the scheme, the state |*Fj.) is 
transformed into I^qds = [c m ic m2 (at™2 +fi^)\mm) + c gl c m2 pf^ x \gm) + c ml c g2 at^ di \rng)}, 
and by choosing yS to satisfy the matching condition af^ + = 0, we have the final state 

l^)ent.QDs = -^=(c ml c g2 \mg) - c gl c m2 \gm)), (21) 



where the normalisation is given by N = \c m ic g2 \ 2 + \c g ic m2 \ 2 . Equation (21) shows that 
1^ )ent.QDs is the desired ideal entangled state \if/~) as long as c m \c g2 = c g \c m2 is satisfied. 
In other words, we can achieve a high fidelity even when the two QDs are not prepared in 
an equal superposition. An equal superposition at the initialisation stage only guarantees a 
higher efficiency 77, but is not required for achieving a high fidelity. This property of the 
scheme provides greater flexibility in the initialisation of the QDs. 



5. Conclusion 



In this work we have described an entanglement generation scheme between two QDs using 
the plasmonic modes of a metal nanoparticle array. Here, we have taken a particular 
nanoparticle structure consisting of two metal nanoparticle arrays which meet at a beam 
splitter of nanoparticles. Although Ohmic loss (energy dissipation) in metals generally 
induces damping in the supported plasmonic systems and limits the use of plasmonics for 
quantum control and state engineering, we have shown a scheme for entanglement generation 
that provides a robust performance against metal loss and the length of the nanoparticle 
arrays. Such robustness overcomes the length limit of nanoparticle arrays for use in nanoscale 
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photonic quantum networks, thus opening up further possibilities for constructing quantum 
plasmonic devices on scales far below the diffraction limit. In addition, in our investigation 
we showed that the fidelity of the entangled QD states was robust against inhomogeneous 
broadenings of the QDs, implying the scheme works even when the QDs have different 
resonant frequencies. Through our analysis, we demonstrated that the robustness of the 
scheme originates from the characteristics of the transmission and reflection amplitudes of 
the system. The robustness against loss comes about as the fidelity mostly depends on the 
behaviour of the reflection coefficient of each arm and the reflection is less sensitive with 
nanoparticle losses since the energy never enters the array. The robustness against detunings 
of the QDs comes about as the operating frequency to be detected is far off-resonant with the 
QDs transitions. Furthermore, we discussed the effect of detection efficiency, the initialisation 
of the QDs, and the use of coherent states as initial states, which reveals additional versatility 
of the scheme in a realistic plasmonic system. The techniques we have employed to 
describe the nanoparticle system in our work may also be helpful in further theoretical 
and experimental studies of plasmonic nanostructures for quantum-control applications and 
probing nanoscale optical phenomena. We hope that the results presented in this paper will 
encourage the use of metal nanoparticle structures with more complex designs for on-chip 
quantum networking. 
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Appendix A. Dipole induced reflection in a metal nanoparticle array 

Here, we provide an analysis of DIR in a single array of metal nanoparticles with one dipole 
(QD), which corresponds to each arm in figure [T] (a). The physical system is depicted in 
figure |A1[ where the n-th nanoparticle is connected to a drain nanowire. The Hamiltonian for 
the nanoparticle array is given by equation ([T]), yielding the Heisenberg equations 

dci r 

— - = -{iO) X + ^ + - igi2&2 + V^^in + ^l^in,l ~ Ud- g(m) (A.l) 

da; Tj . — ~ 

— = -(iiOj + —)aj - igjj.&j.x - ig jJ+l a j+l + V r ^in,j 0' = 2, 3, • • • , n - 1) (A.2) 
clcifj g^ u * I — * 

— = -{ici) n + — + —)ai - ig n - hn a n -i + V&i4i + V r An,i (A.3) 
at 11 

dcf y 

— = -(in + -)& + Uai& z + U (A.4) 
at 2 

where is the natural frequency of the field oscillation at the i-th nanoparticle, g^ k is the 
coupling strength between the fields of the j-th and k-th nanoparticles, Q is the transition 
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(b) Transmission, Reflection, and absorption 
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(c) Transmission, Reflection, and absorption 
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Figure Al. (a) Setup for DIR. A single array of metal nanoparticles with one dipole (QD) 
coupled to the first nanoparticle and a tapered metal nanowire waveguide on the left-hand side 
(source), which focuses light to the end of its tip and excites a localised surface plasmon on the 
adjacent nanoparticle. The excitation propagates across the array of nanoparticles and exits via 
another tapered metal nanowire waveguide on the right-hand side (drain). All metal regions 
have permittivity e m (aS) and dielectric regions have permittivity 6^, as defined in the text. Panels 
(b) and (c) show transmission, reflection, and nanoparticle absorption when S/g np = and 
8 1 gnp = 0-2, respectively, for n = 1. Solid lines are for DIR (QD in |g)), and dashed lines are 
for when there is no DIR (QD in |m». When Aoj « S, DIR occurs, and both transmission and 
loss from nanoparticle dissipation are nearly zero, while reflection is close to one. 



frequency from \g) to \e) of the dipole, / is the vacuum Rabi frequency of the dipole coupled to 
the adjacent nanoparticle, Y t is the damping rate of the /-th nanoparticle, y is the damping rate 
of the dipole, and g s ( d) is the coupling strength of the source (drain) nanowire to its adjacent 
nanoparticle. The operators a] (a*) represent the creation (annihilation) operators associated 
with a dipole-field excitation at /-th nanoparticle, & = \g){e\ 9 and & z = \e)(e\ - \g){g\. The 
operators s m ( 0Ut) and d m ( 0Ut) represent the input (output) annihilation operators of the fields 
of the nanowires, respectively, which satisfy the boundary conditions a\ - -^=(Sin + Sout) an d 
a n - -j=(dm + d out ). The input (output) field of the /-th nanoparticle 's bath is represented by 

bm (out),i, which satisfy the boundary conditions a t = -^(b- m ± + B ou t,i) for j ' = 1, • • • , n, and f d 
denotes the noise operator for the dipole. 

For the above description to be valid all system parameters are restricted by the 
appropriate approximations mentioned in the main text of this paper. We also assume, for 
simplicity, that = ojo for / = g^ = g np , ^ = T for / = and 

gs = = gin-out- The damping rate of the dipole is assumed to be relatively small compared 
to other rates, so that the noise operators, / d , can be neglected in that they do not significantly 
affect quantum coherence of the system. Also, here we consider either longitudinal or 
transverse polarisation along the array according to the direction of electron-charge-density 
oscillations. We impose the weak-excitation regime, where the dipole is assumed to be 
predominantly in the initial state \g). Physically, the weak-excitation limit is valid when a 
single-exctation pulse has a duration that is much longer than the spontaneous lifetime of the 
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dipole ll43ll . Such an approximation has been commonly adopted in many quantum optics 
calculations [|20ll44ll . In this limit, & z can be replaced by its mean value (fr z ) « -1, so that the 
total coupled Heisenberg equations gives rise to the scattering matrix as 

3l(a>) = r s (w)4 t (w) + t s (u)dl t (a>) + f uM (w) (A.5) 
dl(a)) = t d (oj)sl ut (aj) + r d (w)^ ut (w) + f^ d (a)) 9 

where the noise operators for nanoparticles are j^f uM (&0 = Z;4,bi(w)S^ uti (w) anc * ^Ld(^) = 
Z/^d,bi(^)^o Uti (^)- For a given input field at the source (drain) nanowires, the transmission 
and reflection amplitudes are represented by t s ( d ) and r s (d ), respectively, and the coefficient 
b s (d),bi denotes the transition amplitude to a bath mode at the i-th nanoparticle. Here, 
14 (d)(<^)l 2 + ks (d)(<^)l 2 + Zi l^s (d),bi(<^)l 2 < 1? the left hand side does not reach unity as it is 
reduced by the loss of the dipole at the given frequency w for which the dipole interacts 
with the nanoparticle. We define the transmission, reflection, and absorption energy as 
T s (d)(0)) = ks (d)MI 2 , (d)M = ks (d)(^)l 2 , and ^ (d) (co) = Zj l^s (d),bjMI 2 , respectively. 



Appendix A. 1. Single nanoparticle 

We first consider the case of n = 1 in order to provide a basic understanding of the physical 
mechanism for DIR. Arbitrary n is considered in the next section. For a given input field from 
the source nanowire we have 

, A N gin-out(i +*"(<*- Aw)) 

4(Aw) = (A.6) 

J 2 + (| + - Aw))(g in _ out + ^ - /Aw) 

7 2 + (^ + /(5-Aw))(§-/Aw) 

r s (Aw) = — (A.7) 

J 2 + (I + f(<* - Aw))(g in _ out + f - /Aw) 

_ , A N V<?in-out VTo(l + " Aw)) 

Z? sbl (Aw) = ! , (A.8) 

J 2 + (\ + - Aw))(g in _ out + ^ - /Aw) 

where Aw = w-w and 5 = Q-w denote the detunings of the input field from the nanoparticle 
resonance and the dipole resonance from nanoparticle resonance, respectively. If there is no 
dipole in the proximity of the nanoparticle, i.e. J = 0, in the ideal case (T = 0), the field is 
entirely transmitted on resonance (Aw = w - w = 0), i.e. r s (Aw = 0) = and t s (Aw = 0) = 1. 
Once the loss in the metal nanoparticle is included, i.e. T -$> 0, the field is no longer entirely 



transmitted on resonance, i.e., t s (0) ± 1, as shown in figure [AT] (b) as dashed lines. However, 
if there is a dipole in the proximity of the nanoparticle, i.e. 7^0, then the transmission and 
reflection behave differently from the case of no dipole as we discuss next. 



Dipole and single nanoparticle in tune (5 - 0) 



Consider the case where the dipole is resonant with the nanoparticle (5 = 0), equations ( |A.6| )- 
( |A.8| ) can be rewritten as 

?in-out(| - *Aw) 



's(Aw) = - - - / ^ - : (A.9) 

2 



J 2 + a - iAa))(g in - 0Ut + ^ - iAa>) 
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J 2 + a - iAa))W - iAco) 

r s (Aaj) = 2 - 2 - — (A. 10) 

J 2 + (| - /Aa;)(g in _ out + ^ - iAw) 

Z? sbl (Ao;) = . (A.ll) 

J 2 + (| - iAw)(gin_ 0Ut + f - lAtt) 

In this case, when 27 2 /y » g in _ out +r /2 we have 4(0) « and r s (0) « -1, with a 7r phase shift 
on resonance (Aoj = 0), so that the field is totally reflected, similar to the frequency selective 
perfect mirror described by Fan B51l . This is a quantum interference effect called dipole 
induced reflection (DIR). It is based on the same mechanism as dipole induced transparency 
(DIT) originally introduced by Waks et al. EOL in which the cavity-dipole system is driven by 
an external field and the cavity field destructively interferes with the excited state population 
of the dipole. Here, a metal nanoparticle serves as a bad cavity, where cavity decay rate 
T is larger than the dipole decay rate y. Nevertheless, in this bad cavity regime, dipole 
induced reflection (or transmission) still provides strong dispersive properties (as shown in 
figure [ATI (b) as solid lines), which is one of the interesting merits of the original DIT scheme, 
as pointed out by Waks et al EDI . Our plasmonic system for n = 1 is mathematically 
equivalent to driving a double sided cavity with an incident field. 

In order to understand the origin of the dispersive properties, we use the Purcell factor, 
defined as F p = — (proportional to the ratio of the dipole decay rate into the 
nanoparticle, 7, to the bare dipole decay rate, y). In terms of the Purcell factor, the transition 
amplitudes can be rewritten on resonance as t s (Aa> = 0) = t /(F p + 1), r s (Aaj = 0) = 
-(F p - r )/(F p + 1) and b sM (Aoj = 0) = a /(F p + 1), where t = j^j- 2 , r = -j-^fe, 

and ao = g 8in ~°+^ are the transmission, reflection, and absorption amplitudes for a bare 
nanoparticle in the absence of the dipole. One can see that in order to have DIR a large 
Purcell factor is required, i.e. F p » 1. However, as in DIT, we do not need the full normal 
mode splitting condition / > gi n - ut + T /2, known as the high-Q cavity regime. We can 
achieve a large Purcell factor when y «: gi n - ut + T /2 even for much smaller values of This 
is known as the low-Q cavity regime, where the coupling strength / between the nanoparticle 
and the dipole is less than the nanoparticle decay rate T . 

When the nanoparticle plasmonic excitations and the \g)-\e) transition of the dipole are 
in tune (6 = 0), the overall spectrum of the reflection amplitude is always symmetric with 
respect to Aco = 0, at which the local maximum is located, as shown in figure [AT] (b). Here, 
the dipole dissipation rate has a stronger effect on the reflection (transmission) of an on- 
resonance excitation than the nanoparticle dissipation rate does. When y = 0, a plasmon 
is still completely reflected, even with the presence of the nanoparticle loss T , whereas 
when y^0, the dissipation of the dipole interrupts the quantum interference, so that the 
nanoparticle is marginally excited for a plasmon at oj = Q and the loss via the nanoparticle 
increases. Two minima (maxima) in the reflection (transmission) spectrum are located at 
Aco w ±J, in the limit of a large Purcell factor F p » 1, which corresponds to the Rabi-split 
frequencies oj = Q ± J. We find that the linewidth of the broadest transmission peak is the 
nanoparticle linewidth whereas the linewidth of the narrow dip corresponds to the linewidth 
of the dipole dressed by the nanoparticle excitation. This quantum interference that produces 
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the narrow reflection window with a simultaneously strong dispersion results in a significant 
enhancement of the group delay and the possibility of a slowdown or 'storage' of light, as in 
electromagnetic induced transparency [46] . 



Dipole and single nanoparticle detuned (6 ^ 0) 

Consider the case where the dipole is not resonant with the nanoparticle, i.e. 6^0. When 
the dipole transition frequency Q is detuned slightly away from the nanoparticle frequency 
cl>o, the spectrum of the reflection and transmission amplitudes become asymmetric, which 
is known as a Fano resonance, as discussed by Fan [45 j. The local maximum (minimum) 
in the reflection (transmission) spectrum is located at Aoj = 5 in the limit of a large Purcell 
factor F p » 1, regardless of the detuning between the dipole and the nanoparticle, as shown 
by the solid lines of figure |A1| (b) and (c). On the other hand, when the dipole is very 
far detuned from the nanoparticle resonance frequency, the dipole is essentially decoupled 
from the nanoparticle field, so that the reflection (transmission) properties are determined by 
the nanoparticle only, and the reflection (transmission) spectrum dips down to zero at the 
nanoparticle frequency oj = ojq, as shown by the dashed lines in figure |AT] (b) and (c). This 
feature could be exploited to achieve a fast single-excitation switch: for an incoming photon 
with frequency co = cl>o, the transmission is 1 when the dipole is in tune with the nanoparticle 
(Q = ojq), while the transmission is essentially when the dipole is far-detuned. Thus by 
tuning the transition frequency of the dipole, the single-plasmon transport can be regulated 
and the setup acts as a single-plasmon (or photon in the far field) switch, as pointed out in 

Ref. ma. 



Appendix A.2. Arbitrary number of nanoparticles 

We now increase the number of nanoparticles in the array, while keeping the dipole coupled 
to the first nanoparticle, as shown in figure |A1[ As in the case of n = 1, we have transmission, 
reflection, and absorption amplitudes from the scattering matrix given in equation ( |A.6| ). If 



the dipole is decoupled from nanoparticle, the transmission and the reflection amplitudes 
follow the characteristic of a metal nanoparticle array studied in Ref lfT2l . which have different 
dependencies on gi n - ut- We classify their behaviour into two groups: even and odd numbers of 
nanoparticles in the array. As a representative of even and odd numbers of nanoparticles, we 
consider the case ofn = 2 and n = 3. As shown in lfT2ll . for an odd number of nanoparticles, 
there is always a resonance at the natural frequency cl>o, whereas for an even number of 
nanoparticles, a resonance property at the natural frequency ojo depends on the magnitude 
of gin-out/gnp- Thus, the transmission properties as gm-out/gnp is varied at Aoj = are different 
for n = 2 and n = 3 when the dipole is decoupled from nanoparticle, as shown in figure |A2| (a) 
and (b). Here, the transmission (reflection) has a maximum (minimum) value at gin-out /gn P = 2 
for n = 2, whereas the transmission (reflection) rises (drops) quickly as gin-om/gnp is increased 



for n = 3. In addition, in figure |A2| we show the dependence of the transmission and reflection 
with varying amount of metal loss for n = 2 and n = 3. The dependence on the metal loss 
is more sensitive for the transmission than for the reflection, regardless of the parity of the 
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Figure A2. DIR for an arbitrary number of nanoparticles. The transmission T s (oj) (upper) 
and reflection H^co) (lower) at Aoj = is shown as gin-out /gnp is varied for n = 2 (left) and 
n = 3 (right), as representatives of even and odd number of nanoparticles respectively. If the 
dipole is decoupled from its adjacent nanoparticle (solid lines), their behaviour as gin-out /gnp 
is varied becomes different; the transmission (reflection) has a maximum (minimum) value 
at gin-out /gnp = 2 for n = 2, whereas the transmission (reflection) rises (drops) quickly 
as gin-out /gnp is increased for n = 3. Furthermore the transmission and reflection are quite 
sensitive to To/g np changing. On the other hand, if the dipole with a detuning of S/g np = 
is coupled to its adjacent nanoparticle (dashed lines), DIR can be observed when F p » 1 for 
both n - 2 and n - 3. Furthermore, the transmission and reflection are no longer sensitive 
to To/gnp changing; the transmission is nearly zero regardless of gin-out /gnp, and the reflection 
only slightly drops from 1 as gin-out /gnp increases due to the decay of the dipole, y. As the 
number of nanoparticles is increased further, similar trends to those of increasing the metal 
loss are seen, since the effects of increasing n are equivalent to increasing r /g np for a fixed 
number of nanoparticles, as discussed in the main text. 



number of nanoparticles, which is the main reason for the robustness against the metal loss 
mentioned in the main text. As the metal loss is increased, the transmissions become lower; 
for n = 2 (or any even number) they still have a maximum near gin-out /gnp ~ 2, which is 
reflected in the efficiency as shown in figure [4] (c) of the main text. On the other hand, for 
n = 3 (or any odd number) the transmissions are still flat as gin-out /gnp is increased, which 
is also reflected in the efficiency, as shown in figure [4] (d) of the main text. The reflection 
amplitudes shown in A2 (c) and (d) are highly related to the transition amplitude tf^ x shown 
in figure [5] (a) and (b) of the main text, enabling the scheme to be robust to the metal loss. 
On the other hand, if the nanoparticle is coupled to the dipole whose detuning is zero, 



<Vg np = 0, DIR can be observed at Aco = for n = 2 and n = 3 when F p » 1. In figure A2 



while the transmission is nearly zero (dashed lines in (a) and (b)), the reflection (dashed lines 
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012345 012345 012345 



9\v\-ou\lQnp fifin-out/fl'np 5in-out/9'np 

Figure A3. The transmission, T s (oj), as the nanoparticle detuning frequency Aoj/g np and 
input-output coupling gin-out/gnp are varied for arrays of ft = 2, 3, 4, 5, 6, and 7 nanoparticles. 
Regardless of the parity of the nanoparticle number, the early increase of gin-out /gnp enables 
the off-resonant transfer from the source nanowire to the first nanoparticle, whereas its late 
increase leads to strong coupling as if the first nanoparticle becomes the extended 'tip' of the 
nanowire. Thus the large gin-out /gnp implies that the number of nanoparticles is effectively 
reduced to n - 2. Here, DIR occurs at Aoj = 0, when the dipole detuning is given by 
£/gnp = 0, and the transmission becomes nearly zero regardless of gin-out /gnp and the number 
of nanoparticles. 

in (c) and (d)) is nearly one and decreased only slightly with gi n - ut /gnp, which is related to 
the coefficient t™ 8 Sl in figure [5] (a) and (b) of the main text. In addition, when DIR occurs, 
both the transmission and reflection amplitudes are largely insensitive to nanoparticle losses, 
as expected from the case of n = 1 . In figure |A3[ the spectral profile of the transmission is 
presented as Aaj/g np and gin-out/gnp are varied for n = 2, 3, 4, 5, 6, and 7, in the presence 
of metal loss r /g np = 0.1, as in Ref. lfT2l . Note that on resonance (Aa> = 0) DIR can be 
observed regardless of gin-out/gnp and the number of nanoparticles, where the transmission 
becomes nearly zero for Aaj/g np = 0, with the dipole detuning given as 6/g np = 0. 
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